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Abstract 


The  shape  of  the  particle  flux  decline  in  solar  energetic  particle  (SEP)  events  is  of  particular  importance  in  understanding  the 
propagation  of  energetic  particles  in  the  interplanetary  medium.  Power-law  time  profiles  indicate  the  dominance  of  diffusive  prop¬ 
agation,  whereas  exponential-law  decline  emphasizes  convection  transport  and  adiabatic  deceleration.  Values  obtained  theoretically 
for  the  decay  time  in  the  latter  case  are  reasonably  close  to  the  fitted  slopes  in  nearly  half  of  all  events  when  the  solar  wind  speed 
stays  nearly  constant.  Dependencies  of  characteristic  decay  time  t  and  spectral  index  y  on  environmental  plasma  parameters  are 
considered.  Parts  of  exponential-law  declines  when  solar  wind  speed;  (a)  decreases  with  time,  (b)  is  constant,  and  (c)  increases  with 
time  through  the  interval  are  analyzed  separately.  Both  average  values  and  dispersions  of  size  distributions  of  x  for  these  three 
groups  markedly  differ  in  accordance  with  theoretical  expectations. 

©  2005  COSPAR.  Published  by  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

The  time  profiles  of  energetic  particle  fluxes  after  the 
maximum  of  solar  energetic  particle  (SEP)  events  usu¬ 
ally  exhibit  a  several  days  long  decay  phase  of  exponen¬ 
tial  or  power-law  shape.  The  decay  in  each  event  is  going 
on  in  a  different  way,  depending  on  the  state  of  sur¬ 
rounding  interplanetary  medium  and  location  of  accel¬ 
erated  particle  source  relative  to  the  observer.  Different 
propagation  models  result  in  different  shapes  of  the  par¬ 
ticle  flux  declines.  For  impulsive  particle  injection  on  the 
Sun  the  solution  of  the  diffusion  equation  predicts  a  flux 
decline  described  by  power  law  J(t)  °=  t~vl  exp(-r2/4ict) 
at  times  before  convection  becomes  important  ( r  denotes 
distance,  and  k  is  the  diffusion  coefficient,  cf.  Owens, 
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1979).  Except  for  flares,  CME  associated  shocks  may 
be  the  origin  of  accelerated  particles  as  well  (sometimes 
the  only  ones).  During  the  propagation  of  charged  par¬ 
ticles  from  the  Sun,  regions  of  quasi  trapping  of  particles 
can  arise  which  may  either  be  connected  with  the  shock 
front  or  appear  between  the  shock  front  and  strong 
magnetic  fields  at  the  Sun  (Reames  et  al.,  1996;  Kahler, 
1996).  The  decay  of  the  intensity  of  trapped  particles 
also  obeys  to  power  law  (Reames  et  al.,  1996;  Daibog 
et  al.,  2003).  Statistical  studies,  however,  indicate,  that 
the  decay  profiles  of  proton  fluxes  are  predominantly 
exponential.  On  a  sample  of  ~700  events  of  1-5  MeV 
proton  fluxes  only  about  10%  of  the  declines  were 
power-law,  90%  were  exponential,  while  a  small  fraction 
of  exponentials  had  power-law  tails  above  ~20  MeV 
(Daibog  et  al.,  2003). 

Charged  particles  propagating  in  the  solar  wind  plas¬ 
ma,  permanently  present  in  heliosphere,  are  convected 
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and  adiabatically  decelerated  in  the  diverging  magnetic 
field.  Due  to  convection,  in  the  exponential  term  of 
the  profile,  instead  of  -r2  one  obtains  — (r  —  Vt)2,  where 
V  is  solar  wind  speed.  Under  the  most  general  condi¬ 
tions,  from  the  convective  motion  of  accelerated  parti¬ 
cles  it  follows  that  the  x  time  constant  characterising 
decay  should  decrease  with  increasing  V,  and  on  the 
other  hand,  adiabatic  deceleration  should  cause  x  to 
decrease  due  to  energy  spectra  falling  with  energy,  t 
should  also  decrease  with  the  increase  of  spectral  index 
y.  In  case  convection  and  adiabatic  deceleration  terms 
dominate  over  gradients,  scattering,  and  drift  in  the 
transport  equation,  an  exponential  decay  of  flux  is 
obtained:  J  exp(— tlx),  with  a  decay  time  of 

r  =  "ir/2V(2  +  ay),  (1) 

where  a  =  T  +  2mc2l(T  +  me2)  «  2  at  nonrelativistic 
energies  (Forman,  1970;  Jokipii,  1972).  This  can  be 
interpreted  as  adiabatic  deceleration  enhances  the  effect 
of  convection  resulting  in  an  effective  convection  veloc¬ 
ity  V*  =  2(2  +  a y)V  (Owens,  1979  for  impulsive  injec¬ 
tion).  Lee  (2001)  considers  transport  equation  taking 
into  account  adiabatic  deceleration  only  and  obtains 
the  same  dependence  (assuming  a  =  2).  All  results  dis¬ 
cussed  below  refer  to  1  AU,  that  is,  we  cannot  check 
the  variation  of  t  with  r.  The  comparison  of  proton 
fluxes  at  Ulysses  in  the  profiles  of  SEP  events  between 
about  1.3  and  3  AU  near  solar  maximum  with  those  at 
1  AU  (IMP-8)  by  McKibben  et  al.  (2003)  suggests  that 
t  depends  very  little  on  r  above  30  MeV.  However,  espe¬ 
cially  below  ~5  MeV,  the  IMP  profiles  look  steeper, 
which  could  be  attributed  to  a  latitude  effect,  but  an  in¬ 
crease  of  t  with  r  cannot  be  excluded. 

Strictly  speaking,  the  above  formula  is  correct  only  at 
very  late  phase  of  decay,  however,  if  convection  and  adi¬ 
abatic  deceleration  really  take  place,  one  can  expect 
some  correlation  between  t  and  values  of  V  and  y, 
respectively,  after  particle  peak  intensity.  We  consider 
such  correlations  and  compare  them  with  the  formula 
soberly  understanding  the  limits  of  its  application. 


2.  Observations 

The  analysis  of  flux  variations  was  made  on  a  homo¬ 
geneous  data  set  of  2-25  MeV  protons  obtained  by  the 
CPME  instrument  aboard  IMP  8  between  1974  and 
2001.  Time  intervals  during  the  decay  phase  of  SEP 
events  following  flares,  CMEs  and  interplanetary  shocks 
were  examined.  From  the  whole  database  of  more  than 
600  exponential  declines  (see  Daibog  et  al.,  2003)  only 
those  were  selected  which  were  well  approximated  with 
an  exponential  during  a  period  of  at  least  24  h.  As  indi¬ 
cated  in  earlier  analyses,  we  tried  to  minimize  the  influ¬ 
ence  of  shocks  by  only  using  time  intervals  where  the 
particles  directly  related  to  shock  (energetic  storm  parti¬ 


cles)  could  be  neglected.  Events  significantly  “contami¬ 
nated”  by  shock-accelerated  particles  were  omitted  in 
this  analysis.  The  second  requirement  was  regular  solar 
wind  behaviour:  the  variation  of  the  solar  wind  speed  V 
(as  observed  during  the  exponential  periods)  should  ful¬ 
fil  one  of  the  following  criteria:  (A)  V  monotonously  de¬ 
creases  (d  Vldt  <  0),  (B)  V  is  constant  (d  Vldt  =  0)  within 
5%,  or  (C)  V  monotonously  increases  (d  Vldt  >  0). 
Sometimes  A,  B,  and  C  time  intervals  contained  full 
periods  of  exponential  declines  and  vice  versa  (the  latter 
ones  contained  the  formers),  sometimes  they  overlapped 
partly.  As  a  rule,  we  took  into  account  those  periods 
that  overlapped  for  longer  than  12  h.  According  to  these 
criteria  the  numbers  of  declines  in  the  groups  were:  77 


IMP-8  protons 


DOY  1987 


Fig.  1.  Time  profiles  of  protons  in  three  integral  energy  channels  of 
CPME  for  dK/dr<0  (top),  V=  const,  (center),  and  df7dr>0 
(bottom).  Horizontal  bars:  periods  fitted  with  exponentials.  Lower 
curves:  solar  wind  speed. 
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(A),  49  (B)  and  20  (C),  respectively.  The  three  panels  of 
Fig.  1  demonstrate  patterns  of  these  three  groups  of  so¬ 
lar  wind  behaviour  during  exponential  decays.  The  flux 
profiles  were  then  fitted  with  exponentials  and  the  x  val¬ 
ues  determined  from  least-squares  fits.  Here,  we  only 
used  the  2-4.6  MeV  proton  channel. 


3.  Results  and  discussion 

The  solar  wind  speed  values  spread  over  wide  inter¬ 
vals.  Average  values  of  V  are  490,  417  and  498  km/s, 
for  the  three  groups,  respectively.  Parts  of  exponential- 
law  declines  with  different  behavior  of  V  throughout 
the  declines  were  analyzed  separately.  Size  distributions 
of  values  of  x  for  all  three  groups  are  given  in  Fig.  2. 
While  the  values  of  (V)  for  groups  A  and  C  are  practi¬ 
cally  the  same,  the  average  value  of  x  for  the  C  group 
((t)  ~  16.5  h)  markedly  differ  from  both  A  and  B 
((t)  w  19.1  and  19.6  h,  correspondingly). 

In  spite  of  that  all  available  statistics  for  nearly  three 
full  solar  cycles  have  been  used,  the  analysis  of  (V)  and 
x(y)  dependencies  has  limitations.  The  three  panels  of 
Fig.  3  exhibit  scatter  plots  of  t(  V)  for  all  values  of  y  in 
all  groups,  whereas  Fig.  4  gives  the  same  for  r(y)  for 
all  values  of  V.  All  scatter  plots  are  relatively  broad 
which  may  be  explained  by  a  large  dispersion  of  y  for 
x(V)  as  well  as  by  a  large  dispersion  of  V  for  r(y)  plot. 
Furthermore,  the  different  position  (eastward  or  west¬ 
ward)  of  the  observer  relative  to  the  flare  site  (McCrac¬ 
ken  et  al.,  1971)  also  widens  the  distributions.  Another 
possible  cause  of  the  large  scatter  in  Figs.  3  and  4  is 
the  superposition  of  particles  accelerated  by  interplane¬ 
tary  shocks  (especially  for  low  amplitude  events  in 
groups  A  and  C).  The  intensities  of  low-energy 
(>2  MeV)  protons  are  enhanced,  with  a  phase  (often 
lasting  up  to  three  days)  of  ramping  up  before  the  shock 
However,  some  restricting  dependence  can  be  noted  in  a 
sense  that  the  highest  values  of  r  do  correlate  with  small 
values  of  V  whereas  the  highest  values  of  V  correlate 


t  (hr) 

Fig.  2,  Distributions  of  t  values  for  the  three  groups. 


Solar  wind  speed 


Fig.  3.  Dependence  of  r  on  V  for  decays  with  d  Vldt  <  0  (top), 
V -  const,  (center),  and  dVIdt  >  0  (bottom). 

with  small  values  of  x.  The  same  holds  for  the  x(y) 
dependence.  Definitely,  a  tendency  of  a  reverse  x  depen¬ 
dence  on  both  V  and  y  can  be  observed. 

The  complex  dependence  of  x  on  V  and  y  has  been  ex¬ 
plored  by  calculating  correlations  between  experimental 
values  of  x  (texp)  and  those  calculated  according  to  the 
formula  r,h  =  3/72  V  (2  +  ay).  The  two  panels  of  Fig.  5 
display  the  size  distributions  for  the  three  groups.  Aver¬ 
age  values  the  x^x^  ratios  are  equal  to  1.17,  1.04,  and 
0.92,  correspondingly,  thus  the  sign  of  d  VI it  plays  some 
role  in  determining  x.  In  group  A  (xsliplxth)  >  1  practi¬ 
cally  means  as  if  a  higher  effective  value  of  V  should 
be  considered. 

The  two  panels  of  Fig.  5  suggest  that  about  half  of 
the  total  number  of  declines  match  the  value  obtained 
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Fig.  4.  Dependence  of  x  on  y  for  decays  with  d  Vldt  <  0  (top), 
V  =  const,  (center),  and  d  Vldl  >  0  (bottom). 

using  the  formula  (1)  within  an  accuracy  of  about  25%. 
The  dispersions  of  these  size  distributions  are  0.47,  0.42 
and  0.50,  respectively.  Group  C  displays  the  largest 
departure  from  formula  (1).  In  spite  of  that  { V)  has  sim¬ 
ilar  values  in  A  and  C,  and  <y)  is  nearly  the  same  in  all 
groups  (3.0,  3.1,  and  2.8,  respectively),  the  shape  of 
the  distribution  of  C  is  quite  different  from  that  of  A. 
A  possible  explanation  for  such  a  difference  is  that  the 
connection  point  of  the  observer  scans  to  the  east  and 
the  observed  increase  of  the  solar  wind  speed  means  that 
the  faster  solar  wind  presses  the  slower  one  westward  to¬ 
gether  with  the  frozen  magnetic  field  lines  (in  extreme 
cases  corotating  shocks  may  be  formed  at  near  1  AU). 
In  such  a  case  a  very  dynamical  picture  develops  in  a 
region  of  increasing  V  and  probably  conditions  in  this 


15  | - j - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - r 


(b)  T/T„ 

Fig.  5.  Histogram  of  the  distribution  of  x/xxh  values  for  V  -  const,  (left 
panel),  dVIdt  <  0  (upper  line,  right  panel)  and  dV/dt  >  0  (lower  line). 

region  are  not  quasi-stationary.  Here,  both  regular  con¬ 
vection  and  adiabatic  deceleration  regimes  are  violated, 
whereas  in  groups  A  and  C  the  conditions  can  be  consid¬ 
ered  stable.  The  small  number  of  exponential  declines  in 
group  C  and  the  considerably  smaller  value  of  (t)  can 
also  be  a  consequence  of  unstable  conditions  in  inter¬ 
planetary  space  in  the  case  of  increasing  solar  wind 
speed.  However,  in  spite  of  such  unstable  conditions 
sometimes  we  observe  time  profiles,  which  have  smooth 
exponential  character  (see  Fig.  1(a)). 

4.  Conclusions 

The  analysis  of  the  decay  phases  of  nearly  150  SEP 
events  confirmed  the  theoretical  expectation  that  this 
phase  is  predominantly  determined  by  convection  and 
adiabatic  deceleration  in  the  solar  wind  plasma. 

About  the  half  of  the  total  number  of  exponential  de¬ 
clines  analysed  satisfy  theoretical  conditions  and  the  de¬ 
cay  time  falls  within  25%  of  the  theoretical  values.  The 
following  tendencies  in  the  influence  of  solar  wind  speed, 
V,  on  the  rate  of  decline  of  SEP  events  were  revealed: 

•  The  average  values  of  characteristic  decay  times  (t), 
for  groups  of  events  with  decreasing,  constant  and 
increasing  solar  wind  speed  V  are  19.1,  19.6  and 
16.5  h,  respectively. 
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•  t  decreases  with  the  increase  of  solar  wind  speed  and 
of  the  spectral  index  y  in  agreement  with  theoretical 
considerations. 

The  analysis  performed  indicated  that  neither  in  indi¬ 
vidual  events  nor  in  the  entirety  of  exponential  declines 
solar  wind  speed  is  solely  responsible  in  determining  the 
value  of  characteristic  decay  time. 
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I  stavs  pearly  Constant.  Dependencies  of  characteristic  decay  time  t  and  spectral  index  y  on  environmental  plasma  parameters  are 
considered^ 'parts  of  exponential-law  declines  when  solar  wind  speed:  (a)  decreases  with  time  (b)  is  constant,  and  (c)  mcreases  with 
|  time  through  the  interval  are  analyzed  separately.  Both  average  values  and  dispersions  of  size  distributions  of  r  for  these 
I  groups  markedly  differ  in  accordance  with  theoretical  expectations. 
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